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Abstract 

We compute the finite-size corrections to the free energy, internal energy and 
specific heat of the critical two-dimensional spin- 1/2 Ising model on a triangular 
and hexagonal lattices wrapped on a torus. We find the general form of the 
finite-size corrections to these quantities, as well as explicit formulas for the first 
coefficients of each expansion. We analyze the implications of these findings on 
the renormalization-group description of the model. 
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1 Introduction 



It is well-known that phase transitions in statistical-mechanical systems can occur 
only in the infinite- volume limit. In any finite system, all thermodynamic quantities 
(such as the magnetic susceptibility and the specific heat) are analytic functions of 
all parameters (such as the temperature and the magnetic field); but near a critical 
point they display peaks whose height increases and whose width decreases as the 
volume V = L"^ grows, yielding the critical singularities in the limit L — oo. For bulk 
experimental systems (containing V ~ 10^^ particles) the finite-size rounding of the 
phase transition is usually beyond the experimental resolution; but in Monte Carlo 
simulations {V ^ 10^-10^) it is visible and is often the dominant effect. 

Finite-size scaling theory 0) provides a systematic framework for under- 
standing finite-size effects near a critical point. The idea is simple: the only two 
relevant length scales are the system linear size L and the correlation length of 
the bulk system at the same parameters, so everything is controlled by the single ratio 
^oo/^-Q If > Coo, then finite -size effects are negligible; for L ~ ^oo, thermodynamic 
singularities are rounded and obey a scahng Ansatz O ~ L^^Fq^^ooI L) where po is 
a critical exponent and Fq is a scaling function. Finite-size scaling is the basis of the 
powerful phenomenological renormalization group method (see ref. for a review); 
and it is an efficient tool for extrapolating finite-size data coming from Monte Carlo 
simulations so as to obtain accurate results on critical exponents, universal amplitude 
ratios and subleading exponents ^ 0, |^, and references therein] .0 In particular, 
in systems with multiplicative and/or additive logarithmic corrections (as the two- 
dimensional 4-state Potts model |T^), a good understanding of finite-size effects is 
crucial for obtaining reliable estimates of the physically interesting quantities. 

In finite-size-scaling theory for systems with periodic boundary conditions, three 
simplifying assumptions have frequently been made: 

(a) The regular part of the free energy, /j-eg, is independent of the lattice size L |Q 
(except possibly for terms that are exponentially small in L). 

(b) The scaling fields associated to the temperature T and magnetic field h (i.e., Ht 
and /i/i, respectively) are independent of L O . 



(c) The scaling field /i^ associated to the lattice size equals L ^ exactly, with no 
corrections L~^, L~^, ... 0]. 

Moreover, in the nearest-neighbor spin-1/2 two-dimensional Ising model, it was fur- 



ther assumed for many years that there are no irrelevant operators |T^, |T3]; indeed 
this assumption was confirmed numerically through order (T — Tc)^ at least as re- 
gards the bulk behavior of the susceptibility in the isotropic square-lattice Ising model 
||r3|. However, several authors have recently found overwhelming evidence that there 



^ This is true only for systems below the upper critical dimension dc- For Ising models with 
short-range interaction, dc = 4. 

^ Finite-size scaling has also been successfully applied to data coming from transfer-matrix com- 
putations H- 
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are indeed irrelevant operators playing a role in the two-dimensional Ising model 

In particular, for the square-lattice Ising model they have 



g, PI, 0,0,0, 111, H 



found by studying the bulk magnetic susceptibility that there is one irrelevant oper- 
ator contributing to order (T — T^Y and there is (at least) one irrelevant operator 
contributing to order (T — Tc)^. 

An interesting way to test assumptions (a)-(c) and see the effect of the irrelevant 
operators is to compute the asymptotic expansion (in powers of L~^) of the free 
energy and its derivatives with respect to the temperature at the critical point. The 
square-lattice Ising model is the best understood case. 



In a classic paper, Ferdinand and Fisher |2T|] considered the energy and the specific 
heat of the square-lattice Ising model on a torus of width L and aspect ratio p, and 
obtained the first two (resp. three) terms of the large-L asymptotic expansion of 
the energy (resp. specific heat) at fixed x = L(T — T^) [this is the finite-size-scaling 
regime] and fixed p. In particular, at criticality (T = T^) they computed the finite-size 
corrections to both quantities to order L~^. Their results have been improved at the 
critical point by several authors P2, 



Their results can be summarized as 



follows: 



fl\L,p) 

e:^{l,p) 

C2ciL,p) 



/■sq , /2m (P) 

/ bulk "T 2^ T.2m 



I2n 

m=l 

-Eol^i (p) 



m=0 



(1.1a) 

;i.ib) 

(1.1c) 



m=l 



where fc, Ec and Ch,c are respectively the critical free energy, internal energy and 
specific heat.0 

The first important observation is that there are no logarithmic corrections except 
for the specific-heat leading term C^glogL. Secondly, the finite-size corrections are 
integer powers of L~^, which is consistent with irrelevant operators taking integer 
exponents. Furthermore, not all the powers of occur: in the large-L expansion 
of the free energy (resp. internal energy) only even (resp. odd) powers of can 
occur. In the specific-heat expansion all powers of can appear. In addition, the 



■^Janke and Kenna |Q have studied similar expansions for the square-lattice Ising model with 
Branscamp-Kunz boundary conditions. The analytic structure is similar to (1.1) but additional 



terms arise due to the boundary conditions. For instance, there is a term ~ log L/L in the specific 
heat. On the other hand, Lu and Wu |^ studied the critical free energy for the square-lattice Ising 
model on non-orientable surfaces (namely, the Mobius strip and the Klein bottle). They found the 
first terms of the large-L expansion of fdL, p)\ although they did not give details about the analytic 
structure of such expansion. In particular, there is an additional term ~ in the expansion for 
the Mobius strip (due to "surface" effects) which is absent in the Klein bottle. They also explicitly 
showed that the coefficient f2^{p) depends on the boundary conditions [even if the expansion ( l.la| ) 
holds true] . 
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coefficients Cf^ and Ef^ satisfy the relation 

EM _ / -I/V2 for m odd 



Cm{p) [0 for m even 



(1.2) 



The authors of ref. [^] classified (using conformal field theory) all possible irrel- 



evant operators that may occur in the two-dimensional Ising model and found that 
all their results (in the thermodynamic limit and at criticality on a finite torus) can 
be explained in terms of the following conjecture 



Conjecture 1.1 [20, Conjecture (d2)] The only irrelevant operators which ap- 
pear in the two-dimensional nearest-neighbor Ising model are those due to the lattice 
breaking of the rotational symmetry 

In particular, for the square-lattice Ising model the first operator that breaks 
rotational invariance is the spin-four operator + (where here T is the energy- 
momentum operator) whose renormalization-group exponent is y = —2. In ref. pO| 
they showed that this operator can give rise to all the observed corrections in ( |l.l[ ).p| 

In this paper we extend the above results to the triangular and hexagonal lattices. 
We will obtain the large-L asymptotic expansions for the critical free energy, internal 
energy and specific heat for such lattices wrapped on a torus of width L and fixed 
aspect ratio p. The interest of this computation is three- fold. First, we can make 
a new test of Conjecture |1 . 1| . In the triangular lattice, the first irrelevant operator 
(belonging to the identity family) that breaks rotational invariance is -\- with 
y = —6 |2^. If Conjecture ^TT]is true, then several coefficients in the finite-size-scaling 
expansions (|1 . 1|) should vanish. Second, we can directly check whether the ratio (|L 



is universal or not, that is, if ( |1.2| ) depends or not on the microscopic details of the 
lattice. Finally, the asymptotic expansions could be useful to check Monte Carlo 
simulations. 

A ffist study of the triangular-lattice Ising model partition function on a finite 
torus was done by Nash and O'Connor [B^. They obtained (among other interest- 



ing results) the exact expression of such partition function with anisotropic nearest- 
neighbor couplings and extracted its scaling limit. They computed the bulk contri- 
bution to the free energy /bulk and the ffist finite-size correction /2(p)- Here we will 
extend their results at the critical point. 

The main results of this paper can be summarized as follows: 

ML,p) = /buik + f:%^ (1.3a) 

m=l 

EAL,p) = E, + ±^^ (1.3b) 

m=0 



^ A similar finite-size scaUng analysis was carried out for the one-dimensional Ising quantum chain 
which belongs to the same universality class of the two-dimensional Ising model |2^, 29, 31 . 
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oo 



Crn{p) 



CoologL + Co(p) + 



(1.3c) 



m=l 



oo 



Am{,p) 



Ai{p)L + AoologL + Ao{p) + J2 



(1.3d) 



m=l 



f^LiL,P) 



SoologL 



(1.3e) 



where fc (resp. /dog) the third derivative (resp. the logarithmic contribution to 
the fourth derivative) of the free energy with respect to the inverse temperature (3 
evaluated at the critical point. We have also found explicit formulas for the coefficients 



/2(p), hip), Eiip), Er,{p), Coo, C'o(p), Clip), Ciip), C^ip), Aiip), A^o, Ao{p), A(p), 
and Boo (Indeed, = = = Ej = C2 = C3 = A2 = 0). 



Our results on the general analytic structure of the finite-size corrections to these 
models are: 

• The analytic structure of the finite-sizc-scaling corrections of the quantities 
considered here is exactly the same for the triangular and the hexagonal lattices. 

• The finite-size corrections to the free energy, internal energy and specific heat 
are always integer powers of L~^, unmodified by logarithms (except of course for 
the leading logL term in the specific heat). 

• In the finite-size expansion of the free energy, only even integer powers of L^^ 
occur. The only exceptions are the powers L"'' and whose coefficients 



• In the finite-size expansion of the energy, we only find odd integer powers of 
L^^. In this case, the coefficients associated to the powers L^'^ and L^' vanish. 

• In the finite-size expansion of the specific heat, any integer powers of can 
occur, except the terms L^^ and L~^. In addition, the non-zero coefficients of 
the odd powers of in this expansion are proportional to the corresponding 
coefficients in the internal energy expansion as in the square lattice. 

(3) 

• In the finite-size expansion of fc we find that the expected leading term L log L 
is missing, and the actual leading term is simply L. We find that all powers of 
L"^ appear in such expansion, except L~^. 

• In the finite-size expansion of the fourth derivative of the free energy fc^^ we 
find that there is only a logarithmic term ~ log L, even though we expect two 
additional logarithmic contributions of order L log L and log L respectively. 

The above results are very useful to gain new insights on the renormalization- 
group description of the two-dimensional Ising model. Our conclusions on this topic 
are 

• Some irrelevant operators should vanish at criticality. This happens, in partic- 
ular to the less irrelevant one TT with renormalization-group exponent y — —2. 



vanish. 
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In order to give account of all the finite-size corrections, we should include at 



least two irrelevant operators, in agreement with the results of |ll7| , \T8 



The scaling function W{x) (which is the responsible for the logarithmic correc- 
tions to the derivatives of the free energy) vanishes at criticality x = 0. Its first 
derivatives at criticality satisfy 



d^'Wix) 



dx^ 



x=0 



forn = 1,3,4 

l/(A7rv^) forn = 2 



where A = 1 (resp. 2) for the triangular (resp. hexagonal) lattice. These equa- 
tions motivate the conjecture that W{x) = x^/(2A7r\/3). 

The non-linear scaling field associated to the temperature can be computed for 
both lattices and it is given by 

^^t{T) = r-^r' + 0{t') (1.5) 
This result provides a cross-check of the analysis of infinite-volume quantities 



The plan of this paper is as follows: In Section ^ we present our definitions and 
notation. In Sections |^, ||, ^ and ^ we present the computation of the asymptotic 
expansions for the free energy, internal energy, specific heat, and higher derivatives 
of the free energy, respectively. In Section |^ we discuss the consequences of our 
results on the renormalization-group description of the models. In particular, we will 
focus on the irrelevant operators of the model and on the finite-size-scaling functions. 
Finally, in Section |] we present our conclusions and discuss the results. We have 
summarized the technical details in the appendixes: in Appendix ^ we recall the 
Euler-MacLaurin formula, and in Appendix |B| (resp. Appendix P) we collect the 
definitions and properties of the ^-functions (resp. Kronecker's double series). 



2 Basic definitions 

Let us first consider an Ising model on a triangular lattice wrapped on a torus of 
size iV X M at zero magnetic field. The Hamiltonian is given by 

<i,j> 

The partition function is given by 

{a=±l} 

The dual of such triangular lattice is an hexagonal lattice wrapped on a torus of 
size N X M and containing 2NM sites (i.e., the hexagonal lattice can be viewed 
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as a triangular lattice with a two-point basis). The Hamiltonian and the partition 
function of the Ising model on this lattice are also given by ( p.l|) /( ^l2D . 



If one brushes aside some subtleties about boundary conditions, one can relate 
the partition function ( p.2|) of a triangular-lattice Ising model at coupling [3 to the 
partition function of the Ising model on the dual (i.e., hexagonal) lattice at a "dual" 



coupling (3* [33, 34 



Z'^uil^) = Zf^^in 2^-^^*^ (2 sinh 2/3)3^^^/2 (2.3) 
where /3* is defined by 



tanh/3* = e-^'' (2.4) 



Using eq. (|2.3|) and the star-triangle equation we can obtain the critical values of 
the couplings for both models 

{7 log 3 triangular 

1 r- (2.5) 

I log (2 + v3) hexagonal 

However, this argument is strictly valid only in the infinite- volume limit; it gives the 
correct relation 

/*"(/?) = 2f^\n _ 2 log 2 + ^ log(2 sinh 2/3) (2.6) 



between infinite- volume free energies and the correct critical points ( p.5|) , but the 
identity ( p.3|) for finite-lattice partition functions does not in general hold. This 
is because a periodic lattice is non-planar, so that the correct duality formula also 
involves a pair of "homological" modes arising from the two directions of winding 



around the torus [0. Or put it another way: high-temperature graphs that wind 
around the lattice do not necessarily correspond to low-temperature graphs on the 
dual lattice. Therefore, on a finite lattice — which is the subject of this paper — we 
need to be more careful.Q 

We begin by computing the exact partition function of both models on a torus of 



size N X M. One way to do this is by relating the Ising model to a dimer model [^71 



The same computation leading to the square-lattice partition function can be used 



to obtain the hexagonal-lattice partition function by changing the weights of the 
different dimer configurations. Though the triangular-lattice Ising partition function 
cannot be derived from the hexagonal-lattice partition function using duality (|2.3|) , 
for the reasons given above, we can instead use the star-triangle transformation | p5|| . 
Then, the triangular-lattice partition function Z^^^^ is related to the hexagonal-lattice 
partition function ^2MAf (containing twice as much sites) by the formula 

ZtUf^) = RiP)-'"'' Z^l,^0) (2.7) 

^ We thank Alan Sokal for useful clarifications about this point. 
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where the j3 and R{(3) are given by 



sinh2/3 = (2.8a) 
k[(j) sin2p 

Ri(5f = , , ^ , — (2.8b) 

fi:(/5)2 sinh^ 2/3 ^ ' 



and K (which depends on (3 through the parameter v = tanh/3) is equal to 



(1 - v^)'^ 

niP) = , ^ ^ (2.9) 

4^/(1 + t;3)t;3(i + ^;)3 

After straightforward (but lengthy) algebra we find that the partition function for 
both lattices can be written in a very similar way in the ferromagnetic regime: 

Zv.(/?) = ^(2sin2/5)^/2 (2-10) 

a,/3=0,l/2 

where V is the number of spins in the lattice (e.g., V = NM in the triangular lattice 
and V = 2NM in the hexagonal lattice). The functions Za^p{^) are given by 

2 / vr(?2 + a) \ • 2 / T^K^ + P) 



n=0 m=0 ^ ^ ^ 



M 



+ sin' ( ^-^^ - ^ ) + 2sinhV \ (2.11) 



l) triangular 



where the "mass" term /i is given by 

^^ = \''^ ' ^ (2.12) 

2 sinh (5 hexagonal 

The critical point corresponds to the vanishing of the mass, thus giving 

Remark. The fact that the partition function of both lattices depends on the same 
functions Za^isi^) can be explained by noting that the translational symmetry of both 
lattices is the same (i.e., they have the same underlying Bravais lattice). This issue 
explains why the finite-size expansions are so similar in both lattices. 

The functions Za,f^{fj,) can be expanded in powers of /i. In particular, when 
(a,/3) 7^ (0,0) the functions are even in /x, while Zo,o(/^) is an odd function of fi: 

Z^M = Z^,p{0) + ^Z:^^{0)fi' + --- (a,/5)^(0,0) (2.13) 
Zo,o(/i) = /iZ;^(0) + iz::;^(0)/i3 + --- (2.14) 



This is similar to what happens in the square-lattice Ising model |25 
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We are interested in computing the asymptotic expansions for large N and M 
with fixed aspect ratio (e.g. length to width ratio): 



M 



(2.15) 



of the free energy /(/5; A^, p), internal energy E{(3; N, p) and specific heat Ch{i3] N, p) 
at the critical point (3 = (3c- These quantities are defined as follows 



/(/?;iV,p) = -\ogZv{P) 
E{P;N,p) = -^f(P;N,p) 



Ch{P; N,p) 



d(3'- 



■/(/5; N,p) 



(2.16a) 
(2.16b) 

(2.16c) 



In Section ^ we will also consider higher derivatives of the free energy at criticality 

gk 



fi'\N.p) 



dp' 



:f{P;N,p) 



(2.17) 



with k = 3,4. 



Remark. The definition of the specific heat ( p.l6c|) is somewhat non-standard as it 
does not contain the factor 



3 Finite-size-scaling corrections to the free energy 



Let us start with the the basic quantity Za^p ( p. Ill) and write it in the form 

N-lM-l r / / ^ 

. 2 / 7r(r2 + 



= n n ^ 1 + ^^^^ ( 

n=0 m=0 ^ 



'Tc(n + a)\ f2Tc(m + (3) 7r(n + a)\. , , 
cos I — I COS I — — I > (3.1j 



N 



M 



N 



The product over m in ( |3.1|) can be exactly performed with the help of the following 
identity ||3^ : 



M-l 

n 

m=0 



C — A COS 



2n{m + p) 
M 



M 



1 — z_e 



-2m0\ 



(3.2) 



where ( and A are any two real numbers such that |C/A| > 1 and the quantities z± 
are given by 



1 



(3.3a) 
(3.3b) 
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We can finally write Za^p^jj) as 

N-l / 

= 2^^^/2 -Q / ^Qgj^ 2/i + sin^ 0„+„ 

n=0 ^ 

+ \J [cosh 2/i + Sin^ 0„+a] ^ - COS^ ^n.+a^ 

X J] 1 1 - 2_ (n + a, iV, ^)Mg-2.i/3+Mi^„+. I 



A//2 



(3.4) 



n=0 



where we have used the shorthand notations 



z±{k,N,(x) 



cosh 2/i + sin^ 0^ ± W [cosh 2/i + sin^ 0^ 



- cos^ (pk 



cos (pk 



nk 



(3.5a) 
(3.5b) 



Let us now evaluate the functions Zq,^^(0) for (a,/5) ^ (0,0). We follow here the 
procedure used in ref . p5| , which proved to be very efficient for extracting the large- 
asymptotic expansions of the quantities of interest. We first compute the sum 



M 



N-l 



/i = ^ ^ log 1 + sin^ (f)n+a + sin 0n+a\/3 + sin^ 

71=0 L 



M 



N-l 



y^a;i(0n+a) (3.6) 



n=0 



where 



r / 1 . ^ f^p 

Ui(k) = \og 1 + sin2A; + sinA;V3 + sin^A; = Xk +} —X„ (3.7) 



k=2 



The function uji and all its derivatives are integrable over [0, tt], and in addition, 



c.W(7r)-<(0) 



(0), 



/c = 2, 6, 10, 12, 14,... 



otherwise 
We can now use the Euler-MacLaurin summation formula ( |A.6| ) to obtain 

N-l 



(3.8) 



AT ^ 



1 



n=0 



TT 



A 



uJi{(l)n+a) = - j uji{x)dx - ^i^B2{a) 

TT \2™ 52,„+2(a) 



m=l 



(2m + 2)! 



A 



2m+l 



(3.9) 



The first coefficients A^ are 

A = \/3 ; A3 = A7 = ; A5 



16 
73 



A9 = 1792^3; All 



51200 



(3.10) 
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The final result for /i is 

/I = / uji(x)dx Boiaj — npy — rrAom+i 3.11 

" m=l ^ ' 

Let us now consider the quantity /2 

N-l 

/2 = ^ log 1 1 - + a, AT, o^M^-2nip+Mi<t>„+^ I ^3_^2) 

n=0 

We first note that when n + a = N/2, the factor z_{n + a, N, 0) = 0, so this term 
does not contribute to the sum ( p.l2| ). In the other cases z_{n + a,N,0) does not 
vanish and we can use (3.3b) to write ( |3.12| ) as 

N-l 

= J2' log |l — e-A^l°g^+("+"'^'0)-2'rj/3+A^'i<^n+<:« I (3.13) 
n=0 

where Yl' means that we have taken out the term with n + a = N/2 (if such term 
exists) . 

We now proceed as in ref. pSf: we first write log |1 — e~^\ = Relog(l — e~^) and 
then we expand log(l — e~^) as a power series in e~^: 



-i 1 

= - Re V' _(>~MM{logz+in+a,Nfl)-i^„+^)/2+7ri(3] /g 

It is convenient to write the function log z^{k, N, 0) as 

log2;+(fc, A^, 0) =LJ2{4>k) = uji{4>k) - log cos (pk (3.15) 



where ciJi{k) is the function ( p. 71) . We then split the sum over n into two parts: 
n e [0, LiV/2j-l], andn e [LA^/2j , iV- 1]. By making the substitution n - 1 - n 
in the second sum, we finally obtain 



°o lN/2\-l 

-e 



, , L- / -J- 1 



oo Ar-[Ar/2j-l . 

— Rey^ V' _g-2p(A/[M)2(0n + l-a)-i</>n + c]-«7r/3) (3.16) 



p=l 

We now expand the function u)2{k) as a power series in k 

A 



^ ^ ^ 2m + 1 ! ^ ^ 

m=l ^ ' 
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where the are exactly those of the function ui ( p. 10 ). We obtain an expression of 
the form 



oo H [Af/2J-1 



e 



-2p[7rTop(n+o)+i7r/3] 



m=l ^ ' 



\2m+l 



oo ^ Ar-[Ar/2j-l 



e 



-2p[7rrop(n+l+a) — i7r/3] 



p=0 P 



X 



{-.ppg (^)^'" + 1 - (3.18) 



where Tq is a complex number equal to 

r„ = ^ = ^ = e-/. (3.19) 

The next step consists in expanding the exponentials in powers of A^~^. By fol- 
lowing the procedure introduced in [|^, Appendix B] we obtain 



[n + a) 



X exp{— 2p[7rrop(?2 + a) + ivr/?]} 

oo Ar-[Ar/2j-l r °° / \ 2m A 

I ^'^^\N) 2m + 1!^ 



,2m+l 



n ^0 I VA^/ (2m + 1)! 

X exp{-2p[7rrop(n + 1 + a) - ivr/?]} (3.20) 
where the are certain differential operators. The first ones are 

A3 = A7 = (3.21a) 

A5 = A5 (3.21b) 

63 r) 

Ag = + (3.21c) 

All = All (3.21d) 

We can now extend the sum over to n = 00 as the error is exponentially small. 
On the other hand, the contribution of the term with n+a = N/2 is also exponentially 
small, so we can take out this constraint. Then, after rearranging the sums, we obtain 

00 00 

/2 = ^ log 1 1 - e-'^-^[pMn+a)+iP] I _^ ^ log 1 1 _ g-27r[pTo(n+l-a)-i/3] | 

n=0 n=0 

00 2 A 00 00 

E (f) " R<= E E { (" + a)-«e-*'l-<"-)«'" 

m=l p=l n=0 

+ {n + l- c^)2'»+le-2P-[P-o(n+l-a)+i/3]| ^^ ^2) 
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The desired result can be obtained by plugging in ( |B.13|) /( ]"C.2| ): 



/2 



log 



oo 



B,(a) 



m=l 



(2m + 2)!f^^"+^^") 



ReKSl2(^rop)] (3.23) 



where the elliptic ^-function 6a,i3 and the Dedekind's r^-function are defined in Ap- 
pendix the objects Bp{a) are Bernoulli polynomials defined in Appendix and 



^2m+2 Kronecker's double series defined in Appendix 0. Then, the value of 
Zo_^(0) is given by 

logZ„,;3(0) 



ATM, NM 
log2 + 

2 ^ 277 



uji (t) dt + log 



-Tip 



ReK2"i2(zrop) (3.24) 



^_^\NJ (2m + 2)! 
The free energy at the critical point can be computed directly from ( p. 101) : 

UN, M) = -i log 2 + ^ log(2 sinh 2/3J + :^ log Z^,f,{0) (3.25) 



o,/3 



The result ( |3.24| ) means that the free energy for both lattices can be written as 

/2m (p) 



/e(Ar,p) = /b,ik + 5^ 



m=l 



(3.26) 



Thus, only even powers of A^~^ can occur, and in contrast to what happens in the 
square-lattice, we find some even powers whose coefficient vanishes (e.g., /4 = /s = 0). 
The above result agrees with the formula found by Izmailian and Hu for an Ising 
model on a A^x cx) hexagonal (or triangular) lattice with periodic boundary conditions. 
The first coefficients for the triangular lattice are given by 



rtvi 
/bulk 



f2\p) 

/r(p) 

fe\p) 



-log^ + — / cui(t)rft^ 0.8795853861... 
2 v3 27r 



-log 
P 

ptri/ 



\02 


+ 1 


^^3 


+ 


Oil 


2 







fs'ip) = 
Re 



45^3 



\e,\Kl'' + \92\Kl'-' + \9s\Ki 

I ^2 I + I ^3 1 + I ^4 1 



1 i_ 

2 ' 2 



(3.27a) 

(3.27b) 
(3.27c) 

(3.27d) 



where the 9i are the standard 6'-functions defined in ( p.lO|) and the functions A'g 
are given in terms of 6'-functions in (U^^)- As explained in the Appendix ^ all the 
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functions 9i, t] and K^'f^ are evaluated at 2; = and r = ir^p ( |B.11D . The numerical 
values of these coefficients for several values of p can be found in Table |I|. 
The coefficients of the hexagonal-lattice expansion are found to be 



/bulk 



-\og2VQ + ^ I cji(t)rft ^ 1.0250590964 
2 47r Jo 



^hc 
Jim 



(P) 



2/2m(p) 



(3.28a) 
(3.28b) 



The numerical values of the coefficients /^^'^ and /g'^ can be obtained from Table |l| 
with the help of (3.28b). 

Remarks. 1. The values of the bulk critical free energy (|3.27a|) / (|3.28a|) indeed 
coincide with the values obtained from the well-known results in the thermodynamic 
M when (3 = (3^. 



limit 39, 40 



/bulk(/^) 



/bulk(/5) 




JO 



dxdy 

+ log2 

dxdy 
JO ^ 
+ ^log2 




log[cosh^ 2(3 + sinh^ 2(3 — lo{x, y) sinh 



log[l + cosh^ 2(3 — uj{x, y) sinh^ 



(3.29a) 



(3.29b) 



where 



uj[x, y) = cos X + cos y + cos(x — y) 

2. The limiting values of the coefficients /2 and /e as p - 
be [c.f., (B.12)] 



lim /r(p) 

p—>oo 
p—>oo 



V3 



IT 



24 
3l7r5 



(3.30) 

00 are easily found to 

(3.31a) 
(3.31b) 



60480\/3 

The corresponding limiting values for the hexagonal lattice are one half of the above 
values [c.f., (3.28b)]. 

3. Using the properties of the 6'-functions (B.20) / ([B.21| ) and of the functions 
Kq'^ ( |C.5|) / (CTfi) we can easily check that the terms (3.27b)/( |3.27d ) have the correct 



behavior under the transformation < 

/2(P) 

hip) 



M (p l/p). In particular, 

/2(1/P) 

m/p) 
p' 



(3.32) 
(3.33) 



4. From ( |3.24| ) / ( p.21| ) we see that there is in general a non-zero contribution to 
logZa,/3(0) at any order N~'^™' with m > 4. However, we cannot rule out cancellations 
leading to the vanishing of any of the coefficients /2m(p) with m > 5 in (|3.26|) . Similar 
arguments apply to the other large- expansions in the next sections. 
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1 0.636514168294813 

2 0.340929552077890 

3 0.267452513800776 

4 0.242663080213048 

5 0.233284972993438 

6 0.229516370606439 

7 0.227941884733430 

8 0.227265427814348 

9 0.226968542233183 

10 0.226836041806366 

11 0.226776103731001 

12 0.226748689100100 

13 0.226736034656892 

14 0.226730148221756 

15 0.226727392017273 

16 0.226726094176590 

17 0.226725480047888 

18 0.226725188196890 

19 0.226725048974839 

20 0.226724982337581 
oo 0.226724920529277 



0.084178614254145 
0.052778553027830 
0.069393489802385 
0.079193473629707 
0.084621760086675 
0.087503667615417 
0.088999315024483 
0.089766374485276 
0.090157377673790 
0.090356069258841 
0.090456876341452 
0.090507980191895 
0.090533876582898 
0.090533876582898 
0.090553643011573 
0.090557009779414 
0.090558715190155 
0.090559579041296 
0.090560016609581 
0.090560238251165 
0.090560465757793 



Table 1: Values of the coefficients f2\p) and f^^p) for several values of the torus 
aspect ratio p. 



4 Finite-size-scaling corrections to the internal en- 
ergy 

Now we will deal with the internal energy (2.16b). Using (|2.10|) / (|2.13|) we can 
write the critical internal energy as follows: 



£;,(iV,p) = coth2/?, + i^ 



a,/3 



(4.1) 



The derivative dp/ d[3 can be easily computed from eq. ( |2.12| ). Thus, the only unknown 
object is Zqq(O), which can be written as 



N-l . 

+ sin^ (pn + sin 0n \/ 3 + sin^ 0„ 

n=0 ^ 

7V-1 

X JJ |l-2_(n,A^,0)^V^*'^"| 



M/2 



(4.2) 



n=l 
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write (|4.2| ) as 



By noting that the first product is nothing more than /i ( p.6|) with a = 0, we can 

uji{t) dt 



NM NM 
logZ;.o(0) = log2 + log2M + 



vrpA 

N-l 



27r 

oo 



\2m B 



2m+2 



(0) 



m=l 



(2m + 2)! 
+ ^ log 1 1 - (n, AT, 0)^e^''^" | 



A 



2m+l 



(4.3) 



n=l 



The last sum in ( [4 .31 ) is equal to the definition of /2 ( ^.121 ) with a = 0, except for 
the fact that the sum in (|4.3|) starts at n = 1 rather than at n = 0. We can follow step 
by step the same procedure leading to (|3.22| ): the result coincides with ( p.22|) when 
a = except that the first sum in ( p.22| ) now starts at n = 1. Using ( p.l6|) / (|C.2| ) we 
obtain the final result 



, , , NM ^ , NM 

logZ^,o(0) = ^log2 + log2M+^ 
2 Ztx 



ui{t) dt + 2\og\7]{iTQp)\ 



-Tip 



'_^\NJ (2m + 2)! 



ReK. 



0,0 

2m+2 



(4.4) 



This equation implies that the critical internal energy can be written as a power series 
in N-^: 



m=0 



-E,iN,p) = Eo + J2 
For the triangular lattice we find that 
2 



2m+l 



(P) 



]\[2m+l 



3 ^2 ^3 ^4 




\02 


+ 


O3 


+ 


\0a\ 



Et'\p) 



E^ip) = 





6*2^3^4 




15n/3(|^2 


+ 


03 


+ my 



, , i 



\9o\K, 



Re\i\e: 




(4.5) 

(4.6a) 
(4.6b) 
(4.6c) 

(4.6d) 



where we have used ( |B.15| ) / ( |C4| ). The numerical values of these coefficients can be 
found in Table ^ In the hexagonal-lattice case we obtain 



E, 



he 



E, 



he 

2m+l 



(P) 



2 

-^2m+l(p) 

2^3 



(4.7a) 
(4.7b) 
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The numerical values of the coefficients E^^^ and £'5'^ can be obtained from Table ^ 
by using (4.7b). 

Remarks. 1. The limiting values of the coefficients Ei and i?5 as p — > 00 are easily 
found to be [c.f., (B.12)] 

lim Ei{p) = lim E^ip) = (4.8) 

p— >oo p— >oo 

This formula is valid for the triangular and hexagonal lattices. In particular, we expect 
that all the coefficients -E'2m+i(p) will vanish in the limit p ^ 00 due to the existence 
of the factor |^2^3^4| which vanishes exponentially fast. Thus, on an infinitely long 
torus, the internal energy for any finite width N is equal to the bulk value Eq with 
no finite-size corrections. 

2. Using the properties of the 6'-functions (B.20) / ( |B.21|) and of the functions Kq'^ 



(|C3|)/(|C:6D we can easily check that the coefficients Ei and E5 (4.6b)/( ^) /(4.7b) 
have the correct behavior under the transformation p ^ 1/p. In particular, 

M,) ^ ^ (4.9) 
EM = ^ (4.10) 

P^ 



5 Finite-size-scaling corrections to the specific heat 

The specific heat at criticality is given by the following formula 



C 



H,c 



-2 1 d^p 



sinh'2/3c V dp^ 



/3=A 



1 dp 
di3 



E^«,/5(o) 



■E^:,,(o) 

a,/3 



4,0(0) 



\ 



2-1 



/ J 



The main goal of this section is to compute the ratio 

a,l3 



(5.1) 



(5.2) 



where the sums go over (a, (3) 7^ (0, 0). After some algebra, we can write the derivative 
Z"j^{0) as follows 



„ , , AMN 



IT 



V3 



c(l) , oc{2) , r 



(5.3) 
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p 




Ef\p) 


1 


1.017408797595956 


-0.359705063388737 


2 


0.612513647162813 


-0.178088378079924 


3 


0.345040108164264 


-0.168599461543254 


4 


0.185288835745847 


-0.127979167922216 


5 


0.096804501605795 


-0.086206117890971 


6 


0.049827662298672 


-0.054108487929080 


7 


0.025447703091251 


-0.032506071497113 


8 


0.012944169002509 


-0.018975959727317 


9 


0.006570580061525 


-0.010859541565321 


10 


0.003331786807789 


-0.006125084912961 


11 


0.001688570266906 


-0.003416526641140 


12 


0.000855546533105 


-0.001888941526185 


13 


0.000433419665204 


-0.001036828005370 


14 


0.000219555049642 


-0.000565662233279 


15 


0.000111214898315 


-0.000307012198010 


io 


U.UUUUoooo4o4zUoy 


-U.UUUi00oooo4 / iUo 


17 


0.000028535313425 


-0.000089278100310 


18 


0.000014454016292 


-0.000047882460122 


19 


0.000007321388062 


-0.000025601385603 


20 


0.000003708495908 


-0.000013650380771 


oo 









Table 2: Values of the coefficients -Ej'^'(p) and E^\p) for several values of the torus 
aspect ratio p. 



where the sums are given by 



/q 

= ^ E - 

N-l 

Re 



S. 



(2) 
a,/3 



2N 
~2lT 



sm (pn+a V 3 + sin^ 



sin 0„+aA/3 + sin^ 0„+c. ^ 



(5.4a) 



(5.4b) 



The variables 0„+q, and 2;_ = 2;_(?2 + a, N, 0) are given by (3.5) and 6a,o is the usual 
Kronecker's delta. 

The first step is to compute the sum Sa^ ( |5.4a| ). We will follow a procedure similar 
to the one used in ref. p4[ for the square lattice. Let us define the function 

73 



1 1 

- T + 



(5.5) 



sin kV3 + sin^ k k k - % 
This function and all its derivatives are integrable over the interval [0, vr] so we can 
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apply the Euler-MacLaurin formula ( |A.6|) . The final result is 

Af-l 



.n) -W)^^'-' (^-^^ 

where the coefficients 72m-i come from the expansion of ujsik) in powers of k: 

oo ~ 

^3(^) = y^^k^ (5.7a) 



m! 

m=0 



oo , ™ oo 



m=0 m=l ^ ' 

In general, the coefficient 7™, contains two contributions: one comes from the term 
l/(/c — tt) which gives the (trivial) coefficient and the other contribution 

comes from the first two terms in the l.h.s. of ( ^.51) . We will denote by 7^ this 
latter (non-trivial) contribution. In particular, only the coefficients 72m+i with m = 
1, 2, 3, . . . are non-zero. The first non-vanishing coefficients 7^ are 

8 80 448 

73 = Y^; 75 = -^; 77 = — (5.8) 

On the other hand, the value of the integral in ( |5.6| ) is 

I !\,{t)dt = \og^ (5.9) 
In computing the sums Ylin=ScS''^ ~^ ^^^^ result (See e.g., |^) 

E^-'osA' + O^. + ^-Elf^ (5.10) 

n=l k=l 

(where 7^; ~ 0.5772156649 is the Euler constant) and take into account that a = 
0, 1/2. In the simplest case a = we have 



^^"(-' - +'-v-E(^)"*f|--.5;, 

m=l ^ ' 

2{N-1) n) ^^2m{N- 



2m 



(5.11) 



This expression can be simplified by expanding it in powers of A^"^, and then using 
formulas (|A.9|) / (|A.10|) . A further simplification can be made if we take into account 
(5.7b). The final result for a = is 

(iv) = log iv + + log ^ - f: {^y-"- ^^„„_, (5.12) 

m=2 ^ ^' 
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The value for a = 1/2 can be obtained using similar arguments in addition to 
The final result for Sa^ is 



(iV) = log iV + + log ^ - log 4 5^,0 - Yl ( 



m=2 



N) (2m)! 



72m-l (5.13) 



In the above result only enter the non-trivial Taylor coefficients of the function uj^ 

-(2) 

a,/9 



The second step is to compute the sums S^^\ (5.4b). The procedure is similar to 



the ones already done in Sections ^ and We first write 



M ^ g-Aflogz_(n+a,Ar,0) ^ ^-MuJ2{K+a) 



(5.14) 



where the function 002 has been defined in (|3.15|) . Then we split the sum Ylin=Sa 
into two parts: n G [5q,,0) L^/2j — 1] and n G [[A^/2j,iV — 1]. In the second sum we 
perform the change n N — 1 — n and using the properties of cus we arrive at 



0(2) 



TT 



2N 



Re 



-2[M{u)2{4>n+a)-i(t>'n.+a) /2+mj3\ 



L7V/2J-1 

sin0„+.V3 + sinV„+. 1 - e-2[^^(-.(*n..)-^<A„+.)/2+../3] 

Af-L^/2J-1 

/ r^i — i I — r^/ \ 

(5.15) 



+ E 



n=0 



a —>■ 1 — a 



where the second term is the same as the first one with [a, (3) replaced by (1 — a, —(3). 
Now we perform several Taylor expansions: first, we expand the denominator 1 — e"^"^ 
in powers of e"^"^: 



^(2) 



TT 



2N 



Re 



[N/2\ 1 00 _,_2p[Af(a;2(<^n+c)-i</>n+c)/2+7ri/3] 



E E 



sm I 



ha\/3 + 



sm 



Ar_L7V/2j-l 00 g-2p[M(a;2(</'n+i-<.)-i</'n+i-c«)/2-7ri/3] 



+ E E 



^ sm (/) a/3 + sin^0 



(5.16) 



n=0 p= 

Secondly, we expand e~^P'^^'^2/2) g^g ( p,ig|] and finally, we expand the function 



V3 



sin /cV3 + sin^ k 



'^=^^^^ + ^"fc^-^ + ^(2m + l)!^ 

m=l ^ ' 



(5.17) 



in powers of k. After rearranging the series, extending the sums over n to 00 (as the 
error is exponentially small) and using ( [B.14| ) / ( |C.2| ) we obtain 



c(2) 



- Re log Oo^R + 



log2- 



7rp-\/3 



5. 



Q,0 



4E 



TT N 2fc+2 ^2fc+l 



2^ViV; (2A; + 2)! 
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52fc+2(a) - ReK^f+^i^T^^p) 



k,m=l 

oo 



~2~ ^ v a7 



- 



.NJ (2m +1)! 

m=l ^ ' 

where the function W^''^(t) is defined as follows 



(2m + l)!(2fc + 2)! 

2m+l TT/a,/?/- \ 



(5.18) 



:,2m{T{n+a)-l3} 



[n + a) 



n=0 

+ (n + 1 — a 



'I — ^2m{T{n+a)-/3)y 

27ri(r(n+l-a)+/3) 

m ^ 

[\ — g27ri(r{n+l-Q)+/3)-)2 



(5.19) 



The ratio Z^'^(0)/Za^^(0) (|5.3|) can be written as a power series in 



1. 



1 KA^) _ 4 



log + 7s + log 



4^/3 



2 Re log ^, 



a,/3 



4m (P) 
jY2m 



m=2 



(5.20) 



This series contains only even powers of ^ and it starts at ^ (i.e., d 
The first non- vanishing coefficient is 



0). 



(5.21) 



It is worth noticing that the terms with in ( P-3D / ( p.l3| ) / ( |5.18| ) cancel out exactly. 

The computation of the ratio (|5.2| ) is straightforward from ( |3.24|) / (|5.20|) . The 
leading term grows like logA^ and the rest can be expressed as a power series in A^~^ 
where only even powers of A^~^ enter: 



, ^a,/3(0) , oo , . X 

1 ^ ^ AT ^ A ( \ ^ ^ "2m(p) 



(5.22) 



The coefficient associated to A^ ^ vanishes, so the first two non-zero coefficients dm{p) 
are 



do{p) 



_47r3 

"Is" 



, , 4v^ V|e,|Reloge, 
7i5 + log 2- 



TT 



(5.23a) 



+ 



27rp J ( E|g^|Reloggi)(EI^^,/^|Re -^i 
E|ga.;3|Rei^6"'^Relogg^ 

27rp 



1 E 


Reir4"'^) 


73 







E 




El^.. 





(5.23b) 
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where we have denoted by 6i the 6'-functions in the standard notation ( |B.10|) . The 
numerical values of these coefficients can be found in Table |^. 

Remarks. 1. The limiting values of the coefficients c/q and ^4 as p — > oo are easily 
found to be [c.f., (B.12)] 



lim c/o(p) 

p— >oo 



lim c?4(p) 

p^oo 



4 



IT 



V3 



1E + log 



4^3 



TT 



2700\/3 



(5.24a) 
(5.24b) 



2. Using the properties of the 6'-functions (B.20)/( [B.21 ) we can easily verify that 
do{p) has the right behavior under the transformation N ^ M {p ^ 1/p)- 



TC 



V3 



logp + rfo(l/p) = do{p) 



(5.25) 



The behavior of d^lp) under this transformation can be checked numerically to be 
the right one 

= ^ (6.26) 

The specific heat for the triangular and hexagonal lattices can be obtained from 
( ^.1|) and using the results of this section and of Section ^. In particular, we can write 
for both lattices 



ChAN, p) = Coo log N + Co{p) + J2 



m=l 



Cm{p) 



For the triangular lattice the ffist coefficients are given by 

12^3 



Ctri 
00 

C^P) 
Ct\p) 



TT 



<ddo{p)-Q-pEY\pf 
-2EY\p) 

crip) = 

9d,ip)-2pErip)Erip) 
-2Er{p) 



and for the hexagonal lattice the corresponding coefficients are 

2v^ 

2 



Che 
00 



TT 



C^%P) 



■do{p) 



2pEr{p) 



(5.27) 



(5.28a) 

(5.28b) 
(5.28c) 
(5.28d) 
(5.28e) 
(5.28f) 



(5.29a) 
(5.29b) 
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_P do{p) 

1 0.993000152525293 

2 1.205930021583709 

3 1.233520243783654 

4 1.189798214112785 

5 1.134144577781982 

6 1.088416663135744 

7 1.056420958518946 

8 1.035808103247928 

9 1.023167108495542 

10 1.015661512376353 

11 1.011305166086030 

12 1.008818898345350 

13 1.007418256678703 

14 1.006637342854852 

15 1.006205629259680 

16 1.005968648007326 

17 1.005839340365928 

18 1.005769147640350 

19 1.005731215221337 

20 1.005710797002295 
oo 1.005687333437919 



diip) 

-0.034652876469773 
-0.084727027938228 
-0.146295429270869 
-0.167434330275211 
-0.157595538508037 
-0.134596242800881 
-0.110373125092552 
-0.090138741047513 
-0.075083168632513 
-0.064637113563536 
-0.057721974432342 
-0.053296998167206 
-0.050537602364217 
-0.048851578543533 
-0.047838335036578 
-0.047237747222882 
-0.046885886437925 
-0.046681800486309 
-0.046564452961825 
-0.046497492807615 
-0.046411250116879 



Table 3: Values of the coefficients do^p) and di^^p) for several values of the torus 
aspect ratio p. 



Ci^(p) 
C^p) 
C\\p) 



E^ip) 



-d,{p)-ApEr{p)Er{p) 



\ct\p) 



(5.29c) 

(5.29d) 
(5.29e) 

(5.29f) 



The numerical values of the coefficients Cq", C^^^ and Cq'^ can be found in Table ^ 
The values of the coefficients Ci and C5 can be obtained from Table and the value 
can be read from Table § with the help of (|5.29e[ ). 



of Cf 



Remarks. 1. The fact that the coefficients Ci and C3 are proportional respectively 
to El and E^ for the triangular (|5.28| ) and hexagonal ( |5.29|) lattices is not accidental. 
In fact, from ( |5.1| ) / ( |5l2^ we conclude that a// the odd coefficients in the specific-heat 
expansion are proportional to the corresponding coefficients of the internal-energy 
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p 


c'o'Kp) 


c^Kp) 




c^%p) 


1 


1.901880711301990 


0.420058303835062 


0.650313451883665 


2 


4.103024258331998 


-0.326217003543879 


1.017170709722000 


3 


4.744524165326869 


-0.967617404753895 


1.124087360887811 


4 


4.570856116406857 


-1.317204084284658 


1.095142686067810 


5 


4.160445642382107 


-1.334908443794273 


1.026740940397018 


6 


3.780853192640795 


-1.179012991639663 


0.963475532106799 


7 


3.503255527522170 


-0.981777257847270 


0.917209254587028 


8 


3.320932517142029 


-0.807318620952495 


0. 


886822086190338 


9 


3.208115423758777 


-0.674464154921458 


0. 


868019237293130 


10 


3.140842603353847 


-0.581325872529628 


0. 


856807100558975 


11 


3.101715130809265 


-0.519370850894424 


0. 


850285855134878 


12 


3.079361301589703 


-0.479634197647875 


0. 


846560216931617 


13 


3.066761868024451 


-0.454826737355134 


0. 


844460311337408 


14 


3.059735410831789 


-0.439660729459804 


0. 


843289235138631 


15 


3.055850477805811 


-0.430543990999289 


0. 


842641746300969 


io 


o.Uoo III / oiZooo4z 


n /I oc^ 1 on /I o c^n^^o/i o 

-U.4zoioy4zoyooz4o 


0. 


842286296881440 


17 


3.052554049450865 


-0.421972891323654 


0. 


842092341575144 


18 


3.051922325002618 


-0.420136179461406 


0. 


841987054167103 


19 


3.051580935973584 


-0.419080069533795 


0. 


841930155995597 


20 


3.051397172745595 


-0.418477433243636 


0. 


841899528790933 


oo 


3.051186000941275 


-0.417701251051913 


0. 


841864333490213 



Table 4: Values of the coefficients Co"(p), C^\p) and Co'^(p) for several values of 
the torus aspect ratio p. 



expansion. In fact, the proportional constant is given by [See ( ^.1|) /( ^J| )] 



2m+l 



2m+l 



dp 



d'^fi 



d(3^ 



(5.30) 



Indeed, for m = 1, 3 this ratio is indeterminate as both coefficients vanish. 

2. The limiting values of the coefficients Cm{p) as p — oo are easily found to be 
[c.f., (B.12)] 



hm C*"(p) 

p— >oo 

hm cr\p) 

p— >oo 

lim C*"(p) 

p— >oo 



12^3 



TT 



lE + log 



4^/3 



TT 



vr 



2^3 



lim Cf{p) = 

p— +00 

77r3 



300^3 



(5.31a) 
(5.31b) 
(5.31c) 
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lim Cr{p) 

p—>OC 

lim cnp) 

p— >oo 

lim C^'ip) 

p— >oo 



2^3 



TT 



lE + log 



4^3 



7r 



3^3 



lim (p) = 

p— >oo 

77r3 



1800v^ 

3. The behavior of the coefScients Cm{p) under the transformation p 
the expected one 



Co(p) 

Cm{p) 



Coo logp + Co(l/p) 

Cm(l/P) 



for m > 1 



(5.32a) 
(5.32b) 

(5.32c) 
^ 1/p is 

(5.33a) 
(5.33b) 



4. From Table |^ it is clear that C^" should vanish at a value between 1 and 2. 
Actually, due to ( |5.29| e), C^"^ should also vanish at the same value of p. We have 
found numerically that C4 vanishes at 



Pr, 



1.4688897779 



(5.34) 



Indeed, due to the transformation properties of C^i^p) under the transformation p — > 
1/p, C4 also vanishes at p~?^ ~ 0.6807862748. This is similar to what happens in the 
square lattice [E^ . 



6 Higher derivatives of the free energy 



6.1 Finite-size-scaling corrections to /i^^ 

In this section we will consider the third derivative of the free energy (|2.17| ) at 
criticality. Even though this observable is not relevant in practice, its computation is 
interesting as it provides new insights into the finite-size-scaling function W defined 
in Section [^. The observable /i^'* ( ^.17] ) can be written as follows: 

3 



(3) 



8 cosh 2(3c ^ 1 d^P 



sinh=^/3, ' V 



X 



Z^O) 1 (dp\ 



+ 2 



+ 



3 d'^p 



V dp^ 



dp 




13=13, 


f3=l3c 



t.z:a^) / z-o(o) 

E^«,/3(o) VE^«,/3(o) 



^oo(O) 

Y.z.A^) 
2 



(6.1) 



The only unknown object is the derivative Zqq(O), which can be written in the fol- 
lowing way 



z'm 
z'oM 



12MN 



n 



V3 



'^0 



(1) 



90(2) 



(6.2) 
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where the sums were defined in (5.4). By following step by step the procedure 
developed in Section | and leading to (|5.20|) , we can compute the finite-size expansion 
of the ratio (Kl 



A 



1 ^0,0(0) 12 , MMA( 



2m 



(6.3) 



By plugging in ( |6.1|) the above result (|6.3|) and the results already obtained in 
Sections | and ||, we obtain 



/f (iV, p) = A^{p)N + Aoo log N + Aoip) + J2 



m=l 



In this expansion, the coefficient A2 is identically zero. 



<(p) 



1 -16.556352382598901 

2 -16.439945008735128 

3 -6.161165303236093 

4 2.808024778142930 

5 6.829759193109778 

6 7.259238447062151 

7 6.078697042860241 

8 4.522637873985072 

9 3.134987142065343 

10 2.072903487004694 

11 1.324966640033683 

12 0.825436595929018 

13 0.503936252639239 

14 0.302641624812243 

15 0.179285740783859 

16 0.104987265970099 

17 0.060870848500045 

18 0.034988709413814 

19 0.019959520236906 

20 0.011309756764775 
00 



0.588715732188061 
-12.618145549991986 
-16.467144991961212 
-15.425136698441144 
-12.962673854292640 
-10.685119155844771 
-9.019533165133020 
-7.925595102852174 
-7.248692542552664 
-6.845055620123084 
-6.610290784855592 
-6.476167809538218 
-6.400571208146704 
-6.358412464990734 
-6.335102866834869 
-6.322306687731851 
-6.315324296705190 
-6.311533950015708 
-6.309485615841503 
-6.308383036473573 
-6.307116005647652 



(6.4) 



Table 5: Values of the coefficients A^{\p) and v4q"(p) for several values of the torus 
aspect ratio p. 



The most important result contained in ( |6.4| ) is that the coefficient associated to 
the expected leading term ~ L log L vanishes. This is a highly non-trivial fact and 
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1 2.204248900857568 

2 1.260776176148286 

3 0.657422819248581 

4 1.682404023579180 

5 2.442900124268499 

6 2.452406095152028 

7 2.032154004962206 

8 1.508563855204567 

9 1.045164213119409 

10 0.690994951708971 

11 0.441659818805435 

12 0.275146193270710 

13 0.167978851785044 

14 0.100880556816132 

15 0.059761915864342 

16 0.034995755659000 

17 0.020290282882592 

18 0.011662903145113 

19 0.006653173413341 

20 0.003769918921741 
oo 



5.090947597599750 
2.011554369485771 
-0.855894359177661 
-2.151736456053733 
-2.489361833956631 
-2.473169556688759 
-2.378716655165567 
-2.293728702105723 
-2.234638841885926 
-2.197561506291189 
-2.175477057251966 
-2.162714473030565 
-2.155480460188504 
-2.151434956603961 
-2.149195097142150 
-2.147964640036792 
-2.147292993714472 
-2.146928331463238 
-2.146731247829992 
-2.146625156802025 
-2.146503238450814 



Table 6: Values of the coefficients J^^{p) and A^^{p) for several values of the torus 
aspect ratio p. 



we will discuss its physical meaning in Section |^. We have obtained the first four 
non-vanishing coefficients for the triangular lattice 



/itri 



7r 



2pEY\p) 
216 

48 - 54rfo(p) + 6pEj"(p)' 
lQEr\p) 





Re log 9j 







(6.5a) 

(6.5b) 

(6.5c) 
(6.5d) 



where the function (io(p) is defined in (|5.23a| ). The numerical values of J^" and A^q^ 
can be found in Table while the numerical values of A^{^ can be computed from 
Table 0. The coefficients for the hexagonal lattice are 



A\\p) = 2pEf(p) 



he/ 



ApE^%pr + 



36 



71 



V3 





Re log 9j 







— log2|r7 



A he 
^00 



12 

71 



(6.6a) 
(6.6b) 
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= ^-3V3doip)+4V3pE^^{pY (6.6c) 
(p) = 4E^(p) (6.6d) 
and their numerical values can be found in Table |^. 

Remarks 1. The coefficients (|6.5|) / (|6.6|) have the right behavior under the transfor- 
mation p ^ 1/p. In particular, they satisfy 

A(p) = pA,il/p) (6.7a) 
Ao{p) = Aoo\ogp + Aoil/p) (6.7b) 

A.ip) = (6.7c) 
P 

2 - In the limit p — >■ oo, both Ai{p) and Ai{p) go to zero as in this limit 
Ei{p) — > exponentially fast. The limit of the coefficients Aq{p) can be computed 
from (|6jc)/(|gc) and QOlal) 

lim<(p) = 48-^ hi. + log ^ I (6.8a) 

p^co '''' \ TT 

lim<(p) = :^--(lE + \og^] (6.8b) 

3.- From Table ^ we see that ^^^'(p) vanishes at at value pmin between 3 and 4. 
Actually, is zero at 

Pmin,i ~ 3.6249264261 (6.9) 
We also find that v4q"(p) vanishes at 

Pmm,2 ~ 1-0300773853 (6.10) 

In the hexagonal lattice, AY{p) only vanishes in the limit p oo, while Aq'^ is zero 
at 

Pmin,2 ~ 2.6367691963 (6.11) 

6.2 Logarithmic finite-size corrections to /i^^ 

The full finite-size-scaling corrections to the fourth derivative of the free energy 
at criticality /i'^'* are rather cumbersome to compute. However, we can extract with 
much less effort the terms including logarithms. This is what we really need in the 
renormalization-group analysis of Section |^. 

The ffist step is the computation of the full expression for /i^"* in terms of the 
derivatives of the basic objects Za^i^. We should keep only those terms in which 
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Z'^j^{0), Zqq{0) or Z^^^{0) enter. There are three possible contributions 



M) 



ic,log,3 



6 /,,/\2 // ^oo(O) 



/\4 



(6.12a) 
(6.12b) 



E^«./3(0) 
4o(0) 



(6.12c) 



where the derivatives of yU with respect to (3 have been represented for short by 
/i', /i", etc. The first contribution ( |6.12a ) is clearly non-zero and of order logA^. 
The second contribution (6.12b) does not actually contain any logarithm, as the 
logarithmic contributions of Zqq(0)/Zqq(0) and — 3 E ^a,/3(0)'V E ^q,/3(0) cancel out 
exactly [See e.g., (|5.22|) /( |OD ]. The same argument apphes to the second line of 

(iHcD- 

In order to compute the contribution of the first two terms in (|6.12c|) , we have to 
consider the fourth derivative of Za,/^{fi) at /i = when (a,/?) 7^ (0,0). After some 
algebra, we find that the logarithmic contributions to that derivative are 



<Log(0) 



12MA^ 



TT 



V3 



+ —Z^,p{0)\ogN6a,o 



TT 



(6.13) 



where the sums are defined in (5.4). After some more algebra we find that the 
contribution of ( |6.12c| ) does not contain any logarithms. 

In conclusion, we find that the finite-size-scaling expansion for the observable /i^"* 
contains a single logarithmic term 



J c,log 



Boo log AT 



(6.14) 



where Boo can be read from (|6.12a|) . Its numerical value is 



5, 



00 



2736/(7rv^) 
120/(7rV3) 



triangular 
hexagonal 



(6.15) 



The leading term in the large- A^ expansion of /i^'' is expected to be 
also expect a term of order ~ A^. 



A^^, and we 
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7 Irrelevant operators in the two-dimensional Ising 
model 



Let us first collect the main results of the previous sections 



m=l 



EAN,p) = E„ + f %glW (7.1b) 



m=0 



ChAN,p) = CoologiV + Co(p) + f^^^ (7.1c) 



m=l 



/f(iV,p) = A,{p)N + Aoo\ogN + Aoip) + J2^^ (7-ld) 
fci,iN,p) = BoologN (7.1e) 

It is also important to mention that the coefficients f^, fs, E^, E-j, C2, C3, and A2 
vanish. In this Section we will use these results to study the irrelevant^perators 
in the two-dimensional Ising model and the finite-size-scaling function W defined 
below. The results will be applicable to both the triangular and hexagonal lattices 
as the analytic structure of the corresponding asymptotic expansions is the same. 
To our knowledge, there are no predictions based on conformal field theory for the 
hexagonal-lattice Ising model. In this section we will follow basically the notation of 



ref. [20 



Let us start with the basic finite-size-scaling Ansatz for a system defined on a torus 
of linear size L (the aspect ratio is also fixed and plays no role in this discussion), 
zero magnetic field and reduced temperature r [^, Q 

/(r;L) = Mr) + -^W^({/i,(r)L^n) + ^^^({P.M^^^ }) (7-2) 

where /^(t) is a regular function of r and the scaling functions W and W depend 
on the non-linear scaling fields pjir) belonging to the identity and energy conformal 
families. Among them the only relevant field is the one associated to the temperature 
/it(r) (See Table |^. In this Ansatz we have exphcitly used the assumptions (a)-(c) 
introduced in Section |I]. 

The reduced temperature r measures the distance to the critical pointQ and it is 
defined such that t = at P = Pc and r > (resp. r < 0) for /5 < Pc (resp. /? > Pc)- 
In the Ising model on the triangular and hexagonal lattices this parameter takes the 

^ This parameter should not be confused with the torus modular parameter. In this section t 
will mean the reduced temperature, while in the rest of the paper it will denote the usual modular 
parameter. 
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form 



triangular 
hexagonal 



(7.3) 



where as usual v = tanh jS. Under the transformations that map the high-temperature 
phase onto the low-temperature phase and viceversa 

2 

triangular 



1-V 



the reduced temperature simply maps as r 
triangular-lattice case were introduced in 



(7.4) 
hexagonal 

-r. Equations (|7.3|)/(|7^) in the 



The non-linear scaling fields /^^(t) can be written as a power series in r 
/^i(^) = /^i(0) + r^iij + ^rV2,i + • • • 



(7.5) 



and we usually take the normalization fij{0) = 1 for the identity-family fields, and 
fiij = 1 for the energy-family fields (These latter scaling fields are odd under the 
transformation r —>■ — r, thus they satisfy fij{0) = 0). 
As explained in ref. [EO], both scaling functions satisfy 



l^({/i,(-r)(-L)^^}) = W^({/i,(r)L^^}) 



(7.6) 



(and analogously for W). Thus, even (resp. odd) derivatives of W and W with respect 
to r will contain only even (resp. odd) powers of L. This fact explains the structure 
found for the internal-energy and specific-heat expansions: 



ChAL) 



8(3 



dl 

dr 

dl 
dr 



T = 



+ 



T = 



dr 

d(3 



dp 



(7.7a) 
(7.7b) 



r=0 



In particular, (7.7) shows why the odd powers of the specific-heat expansion are 
related to those of the internal energy. We will also make the following assumption, 
which is motivated by the absence of terms L""* log L for any m > in the expansions 

(O 

(d) The scaling function W only depends on the scaling field associated to the 
temperature ^ 

l^({/i,(r)L^^}) = l?(/i,(r)L) (7.8) 

As we know that there are no logarithmic contributions to the free and internal 
energies ( [7.1a]) /(7.1b), the scaling function W should satisfy 



W{0) 



dW{x) 
dx 
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(7.9) 



x=0 



Conformal field theory |20| provides a list of irrelevant operators that may appear 
in the two-dimensional Ising model (See Table |^). By comparing the finite-size-scaling 
Ansatze for the free energy, internal energy and specific heat obtained from ( |7.2| ) to 
the corresponding exact results { ^TT^} we may get new insights about the operator 
content of the model. 



Family 


J 




s y 


[I] 





I 


2 




1 


Tf 


-2 




2 




6 -4 




3 




-6 




4 


Q2Q8 ~^ Q2Q8 


6 -8 




5 


Q12 ~^ Q12 


12 -10 




t 


e 


1 




6 




6 -5 




7 


QaQa 


-7 



Table 7: Operators in the two-dimensional Ising model according to ref. [^. For 
each conformal family, we have listed the primary and quasiprimary fields belonging 
to it. For each scaling field /ij, we show the notation used in ref. pOf, its spin s 
and its renormalizat ion-group exponent y. We have included the most relevant fields 
(i.e., y > —10) with spin s = 6N. We have omitted the conformal family [a] as it is 
irrelevant in this discussion. Only the primary fields I and e are relevant. 



Let us start with the free energy. At the critical point r = this can be written 



as 



(7.10) 



where W depends only on the identity-family fields through the variables Xj = 
/ij(0)L^j, as the energy-family scaling fields vanish at criticality. This expression 
can be Taylor expanded for large L, so we obtain a power series in L^^. The exact 
result ( [r.lap tells us that only corrections of order L^*" can occur, except for the 
terms of order and L~^. From Table |^, we see that the scaling fields TT and /X3 
precisely give corrections of order and to the expansion of ( |7.10| ). Hence, we 
need to impose the conditions 



/iTf (0) = /i3(0) = 

The derivative of the free energy with respect to r can be written as ^0 



dl 
dr 



T=0 



dr 



T=0 



(7.11) 



(7.12a) 



i6H 



dr 



T=0 
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(7.12b) 



Each function Wj{{xk}) can be expanded as we did for the free energy, giving a power 
series in L"^™ with no contribution to orders and L~^. From the exact solution 
(7.1b), we see that only corrections of the type 2,^2m-i ^^^^ appear except for the 
powers and L^^. This implies that the scaling field /ig cannot play any role, thus 

/ii,6 = (7.13) 
The second derivative of the free energy at criticality is given by ||20|| : 



T=0 



T=0 



-21ogL 



dx"^ 



(7.14) 



x=0 



where we have used the standard normalization. The second term in the r.h.s. of 
( [7.14| ) can be written as 



-1 Ly^^''W,,{{xu}) = Wu{{x,}) + ^W„{{x,}) + ■ 



(7.15) 



These two terms alone give all even powers of L ^ except L ^ in agreement with the 
exact expansion ( [7.1c| ). The third term in the r.h.s. of ( [r.l4| ) is equal to 

J2 fi2,Ly^W,i{x,}) = ij;.2,iW^i({xfc}) + j^fi2,2W2i{x,}) + ... (7.16) 



ie[/] 



Again, this a power series containing all even powers of L ^ except L ^ in agreement 
with ([71^ ). 

The coefficient of the leading term should be equal to Coo 



00 



dr 
dp 



W"iO) 



(7.17) 



/3=/3c 



Hence we can determine the numerical value of W"{0) by using (|5.28| a) /( 5.29 a) and 
the definition of r (17.3). The result is 



W"iO) 



l/(27rV3) 



triangular 
hexagonal 



(7.18) 



where we have considered the standard normalization for /it(T). The value ( [7.18|) for 
the triangular lattice agrees with the result obtained in ||20|, eq. (3.34)]. 
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Remarks. 1. The irrelevant scaling fields belonging to the identity family that may 
play a role have non-zero spin (namely, s = 6,12). The spin-zero fields belonging 
to this family should vanish at criticality (e.g., /ij(0) = 0). This result agrees with 
Conjecture |L1[ 

2. The vanishing of the field fii = TT at criticality supports Conjecture (dO) of 
20| . However, our results do not imply their stronger Conjecture (dl): the scaling 

field TT decouples (i.e., fiTTi^) = for all r).[] 

3. In the internal-energy analysis, we concluded that the irrelevant field fi^ should 



vanish at criticality (|7.13|) . This operator has spin six, therefore this result is not 
implied by Conjecture In other words, there are cancellations also in the non- 
scalar sector. On the other hand, we find no constraint on the spin-zero irrelevant 
field yUy. However, if Conjecture |L1| is true, then we should have fiij = 0. 

4. In order to obtain the exact solutions (|7.1| ) we need to include at least two 
irrelevant operators. This result agrees with the findings of [O, [TH| for the square- 
lattice model. It is worth noticing that we can formally obtain the exact solutions 
( [7.1D by including the spin-6 irrelevant scaling field fi2 = T^ + T^ with y = —4 and 
the spin-12 field jUs with y = —10. 



fc{L) 

dl 
dr 



r=0 



T=Q 



^(o) + -^l^(/i2(o)L-^/i5(o)L-l°) 



Oft 
dr 



T=0 

21?" (0) logL + 



+ -Wt{fi2{0)L-\fimL 



-10\ 



+ Wu{fi2{0)L-\fi,{0)L 



T=0 



+ ^/i2,2H^2(/i2(0)L-^/i5(0)L-10) 



/^2,5W^5(/i2(0)L-^/i5(0)L-10) 



(7.19a) 
(7.19b) 



(7.19c) 



Let us now consider the observable /i^"* ( p.l7| ). We are interested here only in the 
terms containing logarithms, which are directly related to derivatives of the scaling 
function W{x). The contribution of this scaling function to this observable can be 
written as 



(3) _ r r ^"'fn\ I 



f^L, = L\ogLW"'{0) 
The exact result ( |6.4| ) shows that 



+ 3\ogLW"iO)^ 



(7.20) 



5W(x) 



(7.21) 



x=0 



^ It is worth mentioning that the authors of pc| | showed by considering the large-distance behavior 
of the triangular-lattice Ising model two-point function that fJ^T fi''') — "(t^)- This result strongly 
supports their Conjecture (dl). 
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This result is consistent with the conjecture put forth by the authors of ref. [2C] who 
claimed that the scaling function W is quadratic in is argument (i.e., W{x) = Ax"^). 

On the other hand, the coefficient of the logarithmic term in /i^"* (i.e., Aqq) is 
proportional to W"{0). This observation provides another way to compute the quan- 
tity W"{0) and a direct mean to test the predictions ( [7.18| ). By using the exact 
results (|6.5| b) / (|6.6|b) / (|7.3|) , we arrive at the same values as in ( [7.18|) supporting the 
correctness of our results. 

Finally, we will discuss the observable fc^\ The contribution of the scaling func- 
tion W to this observable is given by 



^(4) 
c,log 



logL 



(I) 

+31ogL W"{0) 



dr d'^T 
4 h 3 

d(3 8(3^ 



+ 4 



(7.22) 



where all the derivatives of r with respect to P should be evaluated at /3 = /3c 
comparing the above formula to ( |7.1e| ) / ( |6.15| ) , we conclude that 



By 



d^W{x) 



dx"^ 







(7.23) 



x=0 



This result is compatible with W{x) being a quadratic function of x. On the other 
hand, as we know the numerical values of the derivatives of r w.r.t. (3 for the triangular 
and hexagonal lattices, we can use equations ( |6.15| )/ ( 17.221 ) to deduce the value of /is^j. 
The result is the same for both lattices /is^j = —1/4, so the non-linear scaling field /it 
depends on r in the following way 



T 



1 

24' 



+ 0(r5 



(7.24) 



This relation coincides with the function a(r) obtained in ref. pO| for the triangular 
lattice:^ 



a(r) 



1 

r t" + 

24 



47 



161 



(7.25) 



10368 248832 

The equality between a(r) and yUt(T) is important because it provides support to 
Conjecture p..l| : if this conjecture is correct, then both function should coincide pO |. 

We can summarize the results obtained on the scaling function W in the following 
conjecture (which is a natural extension of the conjecture W{x) = x'^/{2tt) for the 
square-lattice model @): 



It is not hard to realize that the function a(T) (|7.25| ) is the same for the hexagonal lattice. The 
key observation is that the free energy for this lattice in the thermodynamic limit (3.29b) can be 
written as 



/he 
bulk 



Jo 



dxdy 
47r2 



log [3 + - uj{x, y)] + cnt. 



where r is given by (O). This equation is equivalent to the definition used in |20| to define a(T) for 
the triangular lattice. 
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Conjecture 7.1 In the Ising model on the trian^lar and hexagonal lattices with 
toroidal boundary conditions, the scaling function W is a function solely of the argu- 
ment X = iJitij)L and this function is equal to 

, . f x^/(27rv^) triangular 

W{x) = l (7.26) 
[ X /(47r\/3) hexagonal 

The coefficient of W should coincide with the constant A obtained in the infinite- 
volume limit analysis of the triangular-lattice model pO|, Eq. (2.34)]. The agreement 



between those coefficients adds support to this conjecture. 



8 Further remarks and conclusions 

We have obtained the asymptotic expansions for the free energy, internal energy, 
specific heat and f^^^ of a critical Ising model on the triangular and hexagonal lattices 
wrapped on a torus of width and aspect ratio p. These expansions are given in 
([7. 1|) . In particular, we have found the exact coefficients /bulk, /4 = /s = 0, /e, 
-^'o. E3 = Ej = 0, E5, Coo, Co, Ci, C2 = C3 = 0, C4, C5, Ai, Aqq, Aq, Ai, A2 = 0, 
and -Boo for both lattices. 

The first important observation is that the analytic structure of the finite-size 
corrections of the observables considered in this paper is the same for the triangular- 
and the hexagonal-lattice models. The reason of this coincidence is that both lattices 
have the same underlying Bravais lattice. This agrees with the physical content of 
Conjecture p..l| : as they have the same rotational symmetry group, they should have 
the same irrelevant operators, hence leading to the same finite-size corrections. 

As it can be seen in (|7.1|) , all the corrections are integer powers of N~^. The 
only exceptions are the logarithmic terms in the specific heat ( |7.1c|) , /i^^ (|7.1d|) , and 
/i^^ ( ra ). In the first case this term is the leading one, while in the other ones 
it is subleading. In the free-energy expansion ( [7. la] ) only even powers of A^~^ can 
occur, while in the internal-energy expansion (7.1b) only odd powers of A^~^ appear. 
In the specific-heat expansion even and odd powers of A^~^ occur. Furthermore, the 
odd coefficients in this latter expansion are proportional to the corresponding odd 
coefficients in the internal-energy expansion. The constant depends on how the mass 
H ( |2.12| ) depends on the temperature ( |5.3CI| ). In the expansion of the observable /^^^ 



we find corrections with all powers of A^~^ except for the term A^~^. Indeed, the 
coefficients fm, Em, Cm and Am do depend on the lattice structure of the model, 
hence they are not universal. 

The fact that E2m+i/C2m+i is a p-independent number for the square lattice (= 
— 1/^/2) p4|, |2^, suggested the idea that this ratio might be universal (e.g., it does 
not depend on the lattice structure) .0 However, our results show that this is not the 

^ Izmailian and Hu (see also computed the finite-size expansion of the free energy 

f{N) ^ /bulk + E^i fk/N^'' and the inverse correlation length C"H^) = E^i h/N^''^^ for a 
critical Ising model on several N x 00 lattices (i.e., square, hexagonal and triangular) with periodic 
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case: 



E. 



2m+l 



(P) 



-1/2 
-73/2 



triangular 
hexagonal 



C2m+1 (p) 

Thus, the proportional constant does depend on the lattice structure, hence it is not 
universal. We can write ( |8.1| ) and the corresponding square-lattice relation ( |1.2| ) in 
an unified way by realizing that the constant is just —1/Eq. 

It is important to noteQ that one key ingredient in this discussion is the fact that 
there is an exact transformation ( [7.4|) mapping the high-temperature phase onto the 
low-temperature phase, so we can define a parameter r ( |7.3| ) transforming as r ^ — r. 
It is not clear whether this transformation exists or not for an Ising model defined 
on a general lattice. However, if that transformation does exist, then we can define 
r so eq. ( |7.6| ) holds, leading to (7.7)/( pTTD . We can summarize all these observations 
in the following conjecture: 

Conjecture 8.1 Let us consider a critical Ising system on a regular two-dimensional 
lattice with toroidal boundary conditions. Let us further assume that there is an exact 
mapping v ^ v' from the high-temperature phase onto the low-temperature phase such 
that T —T. Then, the internal energy and specific heat can be expanded in power 
series of as in f7.1b)/( [7.1d{ ) and the coefficients E^ip) and Cm{p) satisfy 



Em{p) 
CM 



-1/^0 




for m odd 
for m even 



f8.2) 



where Eq is the bulk internal energy [See (7.1b)]. Indeed, we understand that this 
ratio is not defined whenever Em = 0. 



If this conjecture is true, then we could define the expansions 



EciN,p) 
ChAN,p) 



En 



El 



E 

m=0 



E. 



2m+l 



(P) 



]Y2m+l 



Coo logiV + Co(p) + 



m=l 



Cm{p) 



(8.3a) 
(8.3b) 



and then the new ratios would be universal 



Em{p) 

Cm{p) 



-1 





for m odd 
for m even 



•4) 



In Section □ we mentioned that the results contained in this paper could serve also 



to test Monte Carlo simulations. Indeed, the expressions ( p.4|) / ( [4.2|) / ( |5.3|) provide a 
way to compute the exact values of the internal energy and specific heat for any finite 



boundary conditions. They found lattice-dependent coefficients fk and bk, but universal ratios 

;^fc//fc = (22'^--i)/(22^-i-f). 

^°We thank Andrea Pelissetto for useful comments on this matter. 
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torus of size N x M. For very large lattices one could also use the (easier to evaluate) 
asymptotic expansions (7.1b)/ (|7.1c| ). 

On the other hand, by taking the exact values of any observable for fixed aspect 
ratio p and several values of the torus width N, we can check whether the asymptotic 
expansions ( |7. 1| ) are correct or not. In particular, by fitting the exact values to 
the corresponding Ansatz, we can verify whether the numerical coefficients coincide 
with the estimates coming from the fits. We have performed such analysis and we 
have confirmed that the numerical values of the coefficients fm, E^, Cm and Am for 
several values of p coincide with the estimates coming from the fits. In addition, 
this procedure allows us to obtain crude estimates of the next coefficients in each 
expansion. For instance, we obtain for p = 1 (which is the case most frequently 
considered in the literature) the following values 

(8.5a) 
(8.5b) 
(8.5c) 
(8.5d) 

A The Euler-MacLaurin formula 

The Euler-MacLaurin formula is one important tool we need to compute asymp- 
totic series. Here we will use the version of ref. [^, formula 23.1.32]. Let E{x) be a 
function whose first 2n derivatives are continuous in the interval (a, 6). If we divide 
the interval into m equal parts (so that h = (b — a)/m), then we have 

^F{a + kh + ah) = - I F{t)dt 





;i) 


^ 1.932 , 




1) 


^ 0.966 


Eg' 


;i) 


^ -7.821 , 


E^" 


(1) 


^ -2.258 


cr 


;i) 


^ -0.722, 


Cf 


(1) 


^ -0.120 


AT 


;i) 


^ 9.124, 


Af 


:i) 


^ 0.878 



k=l 



— / Bp{a-t)IJ2E^^Ha + kh + th)\dt (A.l) 
^' I k=o J 



where p < 2n, < a < 1, Bn{x) = Bn{x — [_x\) and Bn{x) are the Bernoulli 
polynomials defined in terms of the Bernoulli numbers B^ by 



n 



k=0 



Indeed, -B„(0) = The Bernoulli polynomials satisfy the identity P5| , eq. 23.1.21]: 

5„(l/2) = f-l^ - l) 5„ (A.3) 
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We are mainly interested in sums of the form 

7L-1 



n=0 



27r 



{n + a) 



(A.4) 



The asymptotic expansion of the sum (|A.4|) in the hmit L 00 with 7 fixed can 
be obtained from ( |A.1| ). If we assume that all the derivatives of F{t) are integrable 
over the interval [0, 2717] we can formally extend the sum in (|A.1|) to A; = 00 and 
drop the remainder term [namely, the integral in ( [A. 11 )]. In this case we can write the 
Euler-MacLaurin formula as follows 



n=0 



(n + a, 



1 
2^ 



F{t) dt 



k=i ^ 



Ma) 
k\ 



[F(*^-^)(27r7) -F('=-^)(0)] (A.5) 



In this paper we need the above formula in the particular case L = 2N and 7 = 1/2. 
Then ([A.5|) reads 

^ Af-l 



n=0 



1 /7[_\'' Bkia 



1 r 

- / F{t)dt 



(A.6) 



k=l 



In the computation of the specific heat we also need formula ([A.l| ) in the particular 
case a = and h = 1. In this case we can formally write ( [A.l|) for a function F whose 
derivatives are all integrable over [a, b] in the following form |^ 

En^) = f F{t)dt-\[F{b)-F{a)] 



k=a 



-(2A;) 



where we have used the fact that 

B2k+1 = 



1 

~ 2 




A; = 
A; > 



(A.7) 



As we did in ||2^, we can apply ( [A.7| ) to the function F[x) = with a = and 
6=1. We then obtain the identity 



E 

k=l 



B. 



2k 



2k 



2m 
2k -I 



1 
2 



1 



(A.9) 



If we apply ( |A.7| ) to the case F{x) 
obtain 

m-l 



2m + 1 

^2m-i ■^\^]^ lY^Q same endpoints as before, we 



fc=i 



2m - 1 
2A;- 1 



1 



1 



^11-- 

m 



(A. 10) 
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B Theta functions 

In this appendix we gather all the definitions and properties of the Jacobi's 9- 
functions needed in this paper. We first introduce the object 9ai3{z,T) (a,/? = 
0, 1/2)0 

OaA^, r)=J2 exp L + ^ - a) (z + P - ^)\ (B.l) 

where the nome q is defined in terms of the modular parameter r as follows 

q = e^^^ (B.2) 

Using the identity (proved in |[46||) 



J] [1 + q'^~H] [1 + q^^-H-'] [1 - = ^ g"^ T (B.3) 



n=0 n&L 

we can write (IB. 11) as 



B^^Z.t) = ^(^)gi?2(a)g2.^(l/2-a)(.+/3-l/2) 
oo 

X JJ^ [l _ q2(n+l-a)^2niiz+f3)^ _ ^2(n+a) g-27ri(2+/3) j ^^^^^ 
n=0 

where ri{T) is Dedekind ?7-function 

oo 



71=1 



and B2{a) is the Bernoulli polynomial [c.f., ([A.2 | ) ] 

52(a) = a2-a + i (B.6) 
6 

The relation of the functions 9a,/3 with the usual 6'-functions 9i{z,T) i = 1, . . . , 4 
|7| is the following 

9oAz,t) = ^^(^,r) = -«^(-l)V^*^("+i/2)+"^^("+V2)' (B.7) 

neZ 

9o^i{z,t) = 92{z,t)= ^e2"*^("+i/2)+"^^("+i/2)' (B.8) 

^ig(^,r) = ^4(^,r)= ^(_i)V"^""+"^™' (B.9) 

^1,1 (^,r) = ^3(2, r)= ^e2"'""+"'™' (B.IO) 

This object is almost ident ical to the one introduced in ref. However, this latter one gives 

the wrong sign to 9i{z,t) [c.f. ( B.IO )], although this is not important as we are only interested in 
the case z = where 0i(O, t) — 0. 
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In this paper we will only need these functions evaluated at 2; = and r = irop 
where tq is given by ( p.l9|) . To simplify the notation we will use the shorthands 



Oi = 9i{iTop) = 9i{z = 0,T = iTop) 
V = Vi^'^op) = ?7(r = iTop) 



(B.lla) 
(B.llb) 
(B.llc) 



We also need the limits of the 6'-f unctions in the limit p ^ oo. These limits are given 
by 



lim O^iiT^p) = lim O^^iir^p) = 1 

p— >oo p— »oo 



lim 92{iTop) 



lim 2e-""'"'/^ = 



p— >oo 



p—rCO 



(B.12a) 
(B.12b) 



From eq. (|B.4D we arrive at the following identity valid when {a, (3) ^ (0,0): 



log 



OaA^Top) 



Vii^op) 



7rpRe(ro)52(a) = ^ {log 1 1 - e-^'^l""''^"^^^" 



)~m I 



n=0 



log 1 — e 



- 27r [ro p(n+a) +i/3] 



Another useful relation involving log^a,/3(0, t) is the following 



EE 

n=Sa,o P=l 



n + a 



tTTT 



log^„,/3(r) - ( _ + log2 ) 6a,o 



} (B.13) 



(B.14) 



We have proved this identity by considering each case a,P = 0,1/2 [with 7^ 
(0, 0)] separately and by a careful rearrangement of the corresponding series. 
Dedekind's r^-function satisfies the following identity 



(B.15) 



The analogue of ( |B.13| ) when = (0,0) is given in the particular case r = irop 

by 

00 

^ log |1 - e-2-0pn| ^ log + ^ J^g(^^) (B ig) 

n=l 

We also need the behavior of the 6 functions under the Jacobi transformation 

r^r=-l/r (B.17) 
The result when 2; = is given in ref. HH 



^3(0, r') = {-^r)'/%{0,T) 
^2,4(0, r') = (-try/'e,40,r) 



(B.18a) 
(B.18b) 
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In particular, if r = ir^p where 



< = ^ = - (B.19) 

2 To 



is the complex conjugate of Tq ( |3.19| ), the 6'-functions transforms under ( |B.17D as 
follows 

^3(0, </p) = {Topy/%{0,tTop) (B.20a) 
^2,4(0, zToVp) = (rop)^/'e4,2(0,^rop) (B.20b) 

Finally, we should mention that the absolute value of the above 6'-functions does not 
depend on the sign of ImrQ. Thus, 

\e,{o,tTo/p)\ = \e,{o,ir^/p)\ (B.21) 
C Kronecker's double series 



In this appendix we collect a few properties of the Kronecker's double series 
These series are defined as 



^ ^ m,n£Z ^ ' 

(m,n)^{0,0) 

The basic property we need is the following 



m=l n=0 

+ {n + l- «)2p-le2.im[r(n+l-a)-/3]j ^^ 2) 



in the particular case r = irQp with Tq G C [c.f., ( p.l9| )] and p G M. Eq. ( |C.2|) can 



be easily proved using the same arguments as in ref. [^, Appendix D] where they 
consider the particular case t = ip, p E M.. 

In this paper we also need certain values of the K^'^ obtained in I^S], Appendix 

E] 

„ 1 1 /'7 \ 

(C.3a) 
(C.3b) 
(C.3c) 







^ 1 
30 




- m 


i 


[r) 


= 1 

30 








[r) 


30 




+ ^2'^4 
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K^r) = ^{9l + 9im-6l){9t + et) (C.4a) 



KHr) = ^{01 + et) ^-e^, + etet ] (c.4b) 

•^'0' ^ '04 , Ma / 31 8 4 4 



Krir) = ~^iol + et)i^-et + etet] (c.4c) 
1 1 1 /HI 

i^l'^M = ^ie'2-0t)[-^el-9t9t] (C.4d) 



The behavior of the functions -ft'e under the Jacobi transformation (p. 17] ) can 
be obtained using (B.20) and taking into account that Tq = —1 

K'r{0,tT:;/p) = p'kI'\o,ztop) (C.5a) 
4'°(0,2ro7p) = p'Kly,zTop) (C.5b) 

KpiO,tT^/p) = p'KpiO,tTop) (C.5c) 
Xr(0,</p) = p'K'/{0,zTop) (C.5d) 

Finally we mention that the value of lie Kq'^ does not depend on the sign of Imro: 

ReK^'^iO,tT^/p) = ReKriO,tT,/p) (C.6) 
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